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Federal Research Statement 

This invention was made with government support under Contract No. 
70NAN89H3038 awarded by NIST. The government may have certain rights to the 
invention. 

Background of Invention 

[0001] The present invention relates to a method and system to conduct a 
combinatorial high throughput screening (CHTS) experiment. 

[0002] Combinatorial organic synthesis (COS) is a high throughput screening (HTS) 
method that was developed for pharmaceuticals. COS uses systematic and 
repetitive synthesis to produce diverse molecular entities formed from sets of 
chemical "building blocks." As with traditional research, COS relies on experimental 
synthesis methodology. However instead of synthesizing a single compound, COS 
exploits automation and miniaturization to produce large libraries of compounds 
through successive stages, each of which produces a chemical modification of an 
existing molecule of a preceding stage. Libraries are physical, trackable collections 
of samples resulting from a definable set of the COS process or reaction steps. The 
libraries comprise compounds that can be screened for various activities. 



[0003] 



Combinatorial high throughput screening (CHTS) is an HTS method that 
incorporates characteristics of COS. The CHTS methodology is marked by the 
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search for high order synergies and effects of complex combinations of 
experimental variables through the use of large arrays in which multiple factors 
can be varied through multiple levels. Factors of an experiment can be varied 
within an array (typically formulation variables) and between an array and a 
condition (both formulation and processing variables). Results from the CHTS 
experiment can be used to compare properties of the products in order to discover 
"leads" formulations and/or processing conditions that indicate commercial 
potential. 

[0004] The steps of a CHTS methodology can be broken down into generic operations 
including selecting chemicals to be used in an experiment, introducing the 
chemicals into a formulation system (typically by weighing and dissolving to form 
stock solutions), combining aliquots of the solutions into formulations or mixtures 
in a geometrical array (typically by the use of a pipetting robot), processing the 
array of chemical combinations into products and evaluating the products to 
produce results. 

[0005] Typically, CHTS methodology is characterized by parallel reactions at a micro 
scale. In one aspect, CHTS can be described as a method comprising (A) an 
iteration of steps of (i) selecting a set of reactants, (ii) reacting the set and (iii) 
evaluating a set of products of the reacting step and (B) repeating the iteration of 
steps (i), (ii) and (iii) wherein a successive set of reactants selected for a step (i) is 
chosen as a result of an evaluating step (iii) of a preceding iteration. 

[0006] The study of catalyzed chemical reactions by CHTS involves the investigation of 
a complex experimental space characterized by multiple qualitative and 
quantitative factor levels. Typically, the interactions of a catalyzed chemical 
reaction such as a carbonylation reaction can involve interactions of an order of 6 
or 9 or greater. An investigator must carefully set up a CHTS experiment in order 
to effectively examine such a complex space. Reactant identities and variables, 
process identities and variables and levels of combinations of factors, must be 
chosen to define a space that will provide meaningful results. 

[0007] 

in most instances, an investigator conducts the CHTS experiment for the 
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benefit of a client, who for example, may be a customer from outside the 
investigators company or co-worker from another department within the 
company. In any case, the client attempts to clearly articulate its expectations for 
the experiment to the investigator while at the same time, the investigator 
articulates capabilities and limitations of the CHTS methodology. It is difficult but 
critical to translate the articulations of the client and investigator into an 
experiment definition for the CHTS method. The complexity of a catalyzed 
chemical experimental space makes translation of needs and capabilities into an 
experiment definition even more difficult. There is a need for a method and system 
to conduct an experiment according to specific needs of a client and capabilities of 
the CHTS method. 

Summary of Invention 

[0008] The invention meets this need by a providing a method and system to develop 
an experiment definition for a CHTS experiment. In the method, factors are 
selected for the experiment and interactions among levels of the factors are 
estimated. A probability value of positive interactions is then assigned for each of 
the estimated interactions. A CHTS method is effected on an experimental space 
representing the levels and the probabilities for each interaction are adjusted 
according to results of the CHTS method. 

[0009] The invention also relates to a system for conducting an experiment. The 

system comprises a reactor for effecting a CHTS method on an experimental space 
to produce results and a programmed controller that stores an assigned probability 
value for estimated positive interactions between levels of factors of the 
experimental space and adjusts the probabilities for each interaction according to 
results of the CHTS method. 

Brief Description of Drawings 

[001 0] FIG. 1 is a schematic representation of a system and method for conducting a 
CHTS experiment. 

Detailed Description 
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111] In one embodiment, the invention provides a method and system to permit a 
client and an investigator to confer to develop an experiment definition for a CHTS 
experiment. The method and system can utilize a knowledge matrix as a visual and 
organizational aid to serve as an adjustable definitional model. The matrix model 
can include the factors of the experimental space to be investigated. Determination 
of these factors can require selection of reactant identities and levels and selection 
of process identities and levels and selection of the degrees of combination. For 
example, the experimental factors of the catalyst of a carbonylation reaction can 
be two different metals and a solvent. Levels of one metal may be Fe, Cu, Ni, Pb, 
and Re, of another metal may be V, W, Ce, La and Sn and of the solvent may be 
dimethylformamide (DMFA), dimethylacetamide (DMAA), tetrahydrofuran (THF), 
diglyme (DiGly) or diethylacetamide (DEAA). The model can be set up originally to 
represent an estimation of factor level interactions. The estimation can take the 
form of a probability. The experiment can be conducted and a value of the matrix 
can be adjusted between each iteration of the experiment to represent a 
probability change dictated by the experiment results. 

2] These and other features will become apparent from the drawings and 

following detailed discussion, which byway of example without limitation describe 
preferred embodiments of the present invention. 

3] FIG. 1 is a schematic representation of a system 1 0 and method for conducting 
a CHTS experiment. FIG. 1 shows system 10 including dispensing assembly 12, 
reactor 14, detector 1 6 and controller 1 8. Further shown, is X-Y-Z robotic 
positioning stage 20, which supports array plate 22 with wells 24. The dispensing 
assembly 1 2 includes a battery of pipettes 26 that are controlled by controller 1 8. 
X-Y-Z robotic positioning stage 20 is controlled by controller 1 8 to position wells 
24 of the array plate 22 beneath displacement pipettes 26 for delivery of test 
solutions from reservoirs 28. 



4] 



Controller 1 8 can include a data base repository for storing interaction 
identifications and probability values input by a client or investigator. The 
controller 1 8 also controls aspiration of precursor solution into the battery of 
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pipettes 26 and sequential positioning of the wells 24 of array plate 22 so that a 
prescribed stoichiometry and/or composition of reactant and/or catalyst can be 
delivered to the wells 24. By coordinating activation of the pipettes 26 and 
movement of plate 22 on the robotic X-Y-Z stage 20, a library of materials can be 
generated in a two-dimensional array for use in the CHTS method. Also, the 
controller 1 8 can be used to control sequence of charging of sample to reactor 1 4 
and to control operation of the reactor 14 and the detector 16. Controller 18 can 
be a computer, processor, microprocessor or the like. 

[001 5] An experimental space is defined according to a design that is embodied as a 
program resident in controller 1 8. The design uses input from a client and/or an 
investigator to define interactions and to assign weights that represent 
probabilities that the interactions will be positive. Controller 1 8 translates the 
defined space into a loading specification for array plate 32. Then controller 1 8 
controls the operation of pipettes 26 and stage 20 according to the specification to 
,11 deliver reactant and/or catalyst to the wells 34 of plate 22. 



13 



id 



[001 6] Additionally, the controller 1 8 controls the sequence of charging array plate 22 
into the reactor 14, which is synchronized with operation of detector 1 6. Detector 
1 6 detects products of reaction in the wells 24 of array plate 22 after reaction in 
reactor 14. Detector 16 can utilize chromotography, infra red spectroscopy, mass 
spectroscopy, laser mass spectroscopy, microspectroscopy, NMR or the like to 
determine the constituency of each reaction product. The controller 1 8 uses data 
on the sample charged by the pipettes 26 and on the constituency of reaction 
product for each sample from detector 1 6 to correlate a detected product with at 
least one varying parameter of reaction. 

[0017] 

As an example, if the method and system of FIG.l is applied to study a 
carbonylation catalyst and/or to determine optimum carbonylation reaction 
conditions, the detector 1 6 analyzes the contents of the well for carbonylated 
product. In this case, the detector 1 6 can use Raman spectroscopy. The Raman 
peak is integrated using the analyzer electronics and the resulting data can be 
stored in the controller 1 8. Other analytical methods may be used - for example, 
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Infrared spectrometry, mass spectrometry, headspace gas-liquid chromatography 
and fluorescence detection. 

[001 8] A method of screening complex catalyzed chemical reactions can be conducted 
in the FIG. 1 system 10. According to the method, a client and an investigator 
confer to discuss expectations of the experiment to be conducted in the system 1 0 
and the capability of the system to achieve the expectations. The conference can 
produce a knowledge matrix comprising the experimental space interactions and 
an assigned weighting to each interaction that represent a first estimate of a 
probability that the interaction will be a statistically positive interaction, i.e., that 
the interaction will be a lead. For example, the probabilities can be high, medium 
and low probabilities, represented respectively by numerical weighting values. 
"High, medium and low" mean probabilities that are higher, a medium or lower 
with respect to one another. When three weighting value probabilities are assigned, 
the values can be in respective ranges of about 0.6 to about 0.99 for high, about 
0.2 to about 0.59 for medium and about 0.01 to about 0.1 9 for low. Desirably, the 
respective ranges can be about 0.7 to about 0.9, about 0.2 to about 0.5 and 
about .05 to about 0.1 5. The knowledge matrix is an adjustable definitional model 
that represents the estimated interactions and assigned or adjusted probabilities. 
The model can be a visual organizational aid or the model can be a virtual 
construct resident in a computer database. 

[0019] , . 

Formulations and conditions that represent the interactions are then organized 

according to an experimental design such as a Latin square design or a full 

factorial design. Formulations are prepared according to the design. For example, 

a Latin square design can specify a combination of reactants, catalysts and 

conditions as a multiphase reactant system. In this procedure, a formulation is 

prepared that represents a first reactant system that is at least partially embodied 

in a liquid. Each formulation is loaded as a thin film to a respective well 24 of the 

array plate 22 and the plate 22 is charged into reactor 14. During the subsequent 

reaction, the liquid of the first reactant system embodied is contacted with a 

second reactant system at least partially embodied in a gas. The liquid forms a film 

having a thickness sufficient to allow the reaction rate of the reaction to be 
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essentially independent of the mass transfer rate of the second reactant system 
into the liquid. 

[0020] In one embodiment, the invention is applied to study a process for preparing 
diaryl carbonates. Diaryl carbonates such as diphenyl carbonate can be prepared 
by reaction of hydroxyaromatic compounds such as phenol with oxygen and 
carbon monoxide in the presence of a catalyst composition comprising a Croup 
VIIIB metal such as palladium or a compound thereof, a bromide source such as a 
quaternary ammonium or hexaalkylguanidinium bromide and a polyaniline in 
partially oxidized and partially reduced form. The invention can be applied to 
screen for a catalyst to prepare a diaryl carbonate by carbonylation. 

[002 1 ] Various methods for the preparation of diaryl carbonates by a carbonylation 
reaction of hydroxyaromatic compounds with carbon monoxide and oxygen have 
been disclosed. The carbonylation reaction requires a rather complex catalyst. 
Reference is made, for example, to Chaudhari et al., U.S. Pat. 5,91 7,077. The 
catalyst compositions described therein comprise a Group VIIIB metal (i.e., a metal 
selected from the group consisting of ruthenium, rhodium, palladium, osmium, 
iridium and platinum) or a complex thereof. 

[0022] The catalyst material also includes a bromide source. This may be a quaternary 
ammonium or quaternary phosphonium bromide or a hexaalkylguanidinium 
bromide. The guanidinium salts are often preferred; they include the v , T-bis 
(pentaalkylguanidinium)alkane salts. Salts in which the alkyl groups contain 2-6 
carbon atoms and especially tetra-n-butylammonium bromide and 
hexaethylguanidinium bromide are particularly preferred. 

[0023] 

Other catalytic constituents are necessary in accordance with Chaudhari et al. 
The constituents include inorganic cocatalysts, typically complexes of cobalt(ll) 
salts with organic compounds capable of forming complexes, especially 
pentadentate complexes. Illustrative organic compounds of this type are nitrogen- 
heterocyclic compounds including pyridines, bipyridines, terpyridines, quinolines, 
isoquinolines and biquinolines; aliphatic polyamines such as ethylenediamine and 
tetraalkylethylenediamines; crown ethers; aromatic or aliphatic amine ethers such 



Page 7 of 30 



as cryptanes; and Schiff bases. The especially preferred inorganic cocatalyst in 
many instances is a cobalt(ll) complex with bis-3-(salicylalamino) 
propylmethylamine. 

[0024] Organic cocatalysts may be present. These cocatalysts include various 

terpyridine, phenanthroline, quinoline and isoquinoline compounds including 
2,2':6\2"-terpyridine, 4-methylthio-2,2 , :6\2"-terpyridine and 2,2 , :6 , ,2"-terpyridin 
N-oxide,1,l 0-phenanthroline, 2,4,7,8-tetramethyl-l ,1 0-phenanthroline, 4,7- 
diphenyl-1,10, phenanthroline and 3,4,7,8-tetramethy-1 ,1 0-phenanthroline. The 
terpyridines and especially 2,2':6\2"-terpyridine are preferred. 

[0025] Another catalyst constituent is a polyaniline in partially oxidized and partially 
reduced form. 

[0026] Any hydroxyaromatic compound may be employed. Monohydroxyaromatic 
compounds, such as phenol, the cresols, the xylenols and p-cumylphenol are 
preferred with phenol being most preferred. The method may be employed with 
dihydroxyaromatic compounds such as resorcinol, hydroquinone and 2,2-bis(4- 
hydroxyphenyOpropane or "bisphenol A," whereupon the products are 
polycarbonates. 

[0027] Other reagents in the carbonylation process are oxygen and carbon monoxide, 
which react with the phenol to form the desired diaryl carbonate. 

[0028] The following Example is illustrative and should not be construed as a 

limitation on the scope of the claims unless a limitation is specifically recited. 

[0029] EXAMPLEThis example illustrates an identification of an active and selective 
catalyst for the production of aromatic carbonates. The procedure includes a 
combination of a experimental team weighting procedure and a CHTS method to 
identify a best catalyst from a complex chemical space, where the chemical space 
is defined as an assemblage of possible experimental conditions defined by a set 
of variable parameters such as formulation ingredient identity or amount or 
process parameter such as reaction time, temperature, or pressure. 
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[0030] The chemical space consists of the following TABLE 1 chemical factor levels 
and TABLE 2 processing factor levels: 

TABLE l 





Formulation Type Parameter 
Variation 


Formulation Amount 
Parameter Variation 


Precious metal catalyst 


Held Constant 


Held Constant 


Primary Transition 
Metal Cocatalyst (TM) 


Fe, Cu, Ni, Pb, Re (as their 
acetylacetonates) 


5,10,20,40 (as molar ratios 
to precious metal catalyst) 


Secondary Metal 
Cocatalyst (LM) 


V, W, Ce, La, Sn (as their 
acetylacetonates) 


5,10,20,40 (as molar ratios 
to precious metal catalyst) 


Cosolvent (CS) 


Dimethylformamide (DMFA), 
Dimethylacetamide (DMAA), 
Diethyl acetamide (DEAA), 
Tetrahydrofuran (THF), 
Diglyme (DiGly) 


50,100,200,400 (as molar 
ratios to precious metal 
catalyst) 


Hydroxyaromatic 
compound 


Held constant 


Sufficient added to achieve 
constant sample volume 



[0031] 



Process parameters are shown in TABLE 2: 
TABLE 2 



Process Parameter 


Parameter Variation 


Temperature 


Constant atl00°C 


Pressure 


Constant at 1500 psig 



[0032] Pre-test estimates of interactions among factor levels are postulated at a 
meeting between a customer and investigators. The estimates are assigned 
probability values, which are expressed in the following knowledge matrix TABLE 
3. The probabilities are constrained to three possible values, 0.8, 0.3 and 0.1 f 
which express high, medium, and low probabilities. Probabilities of 0.0 and 1.0 are 
excluded from off-diagonal cells since these probabilities imply complete 
knowledge. The matrix is symmetrical around the main diagonal, since the 
probability of A interacting with B is the same as the probability of B interacting 
with A. 









TABLE 3 










TM 


TM 


LM 


LM 


CS 


CS 


TMType 


Type 


Amount Type 


Amount 


Type 


Amount 


1 


0.8 


0.3 


0.3 


0.3 
0.3 


0.3 
0.1 


TM Amount 


0.8 


1 


0.3 


0.1 


LM Type 


0.3 


0.3 


1 


0.8 


0.3 


0.1 


LM Amount 


03 


0.1 


0.8 


1 


0.1 


0.1 


CS Type 


0.3 


0.3 


0.3 


0.1 


1 


0.8 


CS Amount 


0.3 


0.1 


0.1 


0.1 


0.8 


1 
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[0033] 



The matrix information is loaded into a computer database. The computer 
defines a full factorial experiment according to two factor interactions between 
levels as shown in TABLE 4. The computer also controls a dispensing assembly and 
loading robot to load experimental array trays and a reactor to conduct a CHTS 
experiment. In the experiment, catalyzed mixtures are made up in phenol solvent 
using the concentrations of each component as given in the rows of TABLE 4. The 
total volume of each catalyzed mixture is 1 .0 ml. From each mixture, a 25 
microliter aliquot is dispensed into a 2 ml reaction vial, forming a film on the 
bottom. The vials are grouped in array plates by process conditions (as specified in 
the TABLE 2 Pressure and Temperature columns) and each array plate is loaded 
into a high pressure autoclave and subjected to the reaction conditions specified. 
At the end of the reaction time, the reactor is cooled and depressurized and the 
contents of each vial are analyzed for diphenyl carbonate product using a gas 
chromatographic method. Performance is expressed numerically as a catalyst 
turnover number or TON. TON is defined as the number of moles of aromatic 
carbonate produced per mole of Palladium catalyst charged. This is shown in 
column TON of TABLE 4. 
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TABLE 4 



TMType LMType CSType 


TMAmt 


TMAmt 


V^O /Villi 


TOM 


Ni 


v 


Lsl\Jiy 


10 


< 


9ftft 


1U54 


Re 




HP A A 


10 


1ft 


9ftft 


ijy4 


Ni 




L/IYLrA 


10 


4ft 


AfY\ 
£ fV\) 


1221 


Ni 


Sn 


FIT? A A 


40 


< 

J 


<n 


lo97 


Fe 


T st 


L/MAA 


10 


9ft 


9ftft 


949 


Pb 




I-HtiAA 


4ft 


1ft 


<n 


2317 


Fe 




TTJTJ 


4ft 


9ft 


inn 

IUU 


792 


Cu 


V 


txtc 
IrLr 


1ft 


1ft 


ZUU 


1054 


Cu 


on 




2ft 


1ft 


inn 
IUU 


1 At O 

1058 


Cu 


T a 


LMVLAA 


1ft 


4ft 


^n 
3U 


1081 




v_e 


PlIOTDA 

lJMrA 


K 

J 


4ft 


inn 
IUU 


1058 


Re 

£\.C 


V 

V 


THF 


s 

J 


c 

3 


4UU 


1074 


Cu 


W 

TV 


TWA A A 






Ann 


1125 


Cu 




TUC 

lxir 




1ft 


onn 
2UU 


1111 
1111 


Ni 


T j» 


r\vrc a 


2ft 


K 


<n 

DU 


1358 


Fe 




I/lVlrA 


1ft 


9ft 


<n 


955 


Pb 


V 


r»p A A 


1ft 


c 

J 


inn 
IUU 


1040 


Cu 


v 


nmjA 

UiYirA 


20 


4ft 


3U 


1UV2 


Re 


v 


DMA A 


10 


1ft 


inn 
IUU 


lUoU 


Re 


v 


TWA A 


10 


1ft 


jU 


i n>iQ 
1U49 


Pb 


V 


DEAA 


20 


10 


4ftft 




Pb 


Ce 


THF 


10 


10 


100 


1248 


Fe 


W 


DMAA 


40 


40 


50 


914 


Ni 


La 


DMAA 


5 


20 


50 


1069 


Fe 


La 


DEAA 


5 


20 


400 


1069 


Cu 


Sn 


DiGly 


20 


40 


200 


1114 


Cu 


W 


DMFA 


40 


10 


200 


1105 


Pb 


Sn 


DMAA 


10 


40 


50 


1511 


Fe 


V 


THF 


40 


10 


400 


1067 


Re 


W 


DiGly 


5 


10 


400 


1034 


Cu 


W 


THF 


20 


10 


400 


1041 


Pb 


La 


THF 


10 


5 


50 


1371 


Pb 


V 


DMFA 


20 


10 


100 


1056 


Ni 


W 


THF 


20 


5 


200 


1136 


Pb 


Sn 


DiGly 


10 


10 


200 


1499 


Re 


W 


DMFA 


40 


5 


50 


1535 


Ni 


Ce 


DEAA 


10 


20 


200 


1164 


Re 


Ce 


DEAA 


40 


20 


50 


1959 


Re 


La 


DiGly 


5 


40 


200 


1077 


Pb 


La 


DEAA 


5 


40 


50 


1108 


Re 


Sn 


DiGly 


10 


40 


400 


1660 


Ni 


V 


DMFA 


5 


5 


100 


1083 
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Re 


La 


i'HK 


Re 


Sn 




Fe 


Cc 


DMFA 


Re 


W 




Re 






Ni 


Ce 


DEAA 


Re 


W 


U\\Jiy 


Pb 


W 


THF 


Ni 


Ce 


THF 


Cu 


v 

V 


THF 


Fe 


Sn 


HFAA 


Fe 


La 


DMA A 


Re 


v 

V 


TUT3 


Ph 


W 


DlCjly 


Fe 


W 

VY 




Fe 


on 




Ni 


v 




Ni 


Oil 


T\»yr A A 


Re 


on 


TUTJ 


Ni 


W 


i/ioiy 


Pb 


Sn 

on 


TYB A A 


Cu 




HPT TT? 

Irir 


Ni 


w 


uioiy 


Ni 


Ce 


DFAA 


Re 


v 


uivjiy 


Cu 


v 


fulfil*/ 


Cu 


Sn 


UMrA 


Ni 


Ce 




Pb 


v 




Fe 


W 


DEAA 


Fe 


w 


THF 


Re 


Ce 


DMAA 


Ni 


W 


DMFA 


Cu 


La 


DiGly 


Pb 


Ce 


DMFA 


Pb 


V 


DMAA 


Pb 


V 


DMAA 


Ni 


Sn 


DMAA 


Ni 


V 


THF 


Re 


La 


DMFA 


Pb 


Ce 


DiGly 


Cu 


Sn 


DMAA 


Fe 


Ce 


THF 


Pb 


V 


DMFA 


Re 


Sn 


DMFA 


Pb 


V 


THF 



4U 


5 


200 


2396 


on 


1 A 

10 


100 


2291 


c 

J 


40 


200 


1029 


40 


5 


400 


1538 


1 A 

10 


20 


100 


1417 




10 


400 


1251 


oa 
20 


5 


100 


1376 


5 


20 


50 


1058 


in 
10 


5 


100 


1236 


OA 

20 


40 


100 


1078 


10 


10 


400 


837 


oa 

20 


20 


50 


805 


40 


20 


50 


1076 


10 


20 


100 


1194 


10 


40 


200 


1017 


10 


5 


50 


857 


20 


20 


100 


1065 


40 


40 


400 


1645 


40 


10 


100 


2878 


Af\ 

40 


40 


100 


1173 


5 


5 


400 


1080 


try 
40 


5 


50 


1038 


20 


10 


200 


1215 


OA 

20 


5 


200 


1275 


OA 

20 


5 


50 


1085 


1 A 

10 


20 


400 


1046 


1 A 

10 


5 


400 


1093 


f 

J 


5 


50 


1069 


* 


An 


400 


1039 


40 


5 


400 


936 


10 


10 


100 


1043 


20 


5 


100 


1705 


20 


20 


100 


1187 


5 


10 


50 


1098 


20 


40 


50 


1458 


5 


10 


200 


1113 


40 


20 


50 


1072 


5 


5 


100 


1089 


20 


40 


50 


1092 


10 


20 


200 


1531 


5 


10 


50 


1067 


5 


20 


200 


1034 


5 


40 


400 


1105 


10 


40 


200 


1110 


5 


10 


50 


1078 


5 


20 


200 


1136 
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Ni 


La 


DEAA 


in 


< 


Fe 


Sn 


THF 


«■ 

. J 


5 


Pb 


La 


DMA A 


e 

J 


40 


Cu 


Ce 




on 


5 


Ni 


W 


DP A A 


e 
J 


40 


Pb 




■uioiy 


5 


5 


Pb 


Sn 


TT-TF 


on 


40 


Cu 




r\\jrc A 


10 


10 


Re 


Sn 
on 


uiuiy 


5 


20 


Re 


TV 


l rlr 


10 


40 


Ph 


w 


TM2 A A' 


5 


10 


Cu 


V 

V 


y\rj a A 


10 


5 


Cu 




UrnrJ\ 


40 


20 


rc 


on 


DMrA 


5 


40 


Pb 




JJJVLrA 


40 


20 


Pn 

v-u 


W 


THF 


5 


5 


Fe 


v 

V 


UmrA 


40 


40 


Ni 
lit 


*-e 


THF 


10 


10 


Ph 


v_e 


r\\i a A 


20 


40 


Fe 


Oil 


ntjA a 


20 


5 


Ni 


Sn 

on 


Irlr 


10 


40 


Cu 




ui\Jiy 


10 


40 


Pb 




TWA A A 
JJMAA 


40 


5 


Fe 


v 




5 


10 


Pb 


w 


DMA A 


OA 


10 


Re 


Ce 


DMA A 


5 


40 


Fe 


La 


DMA A 


on 


10 


Re 


V 




5 


40 


Fe 


Sn 


THF 


< 




Ni 


V 


DiGly 


40 


10 


Ni 


V 


DMAA 


20 


5 


Fe 


La 


DiGly 


20 


40 


Cu 


La 


DEAA 


5 


5 


Re 


La 


DiGly 


20 


20 


Pb 


Ce 


DEAA 


40 


5 


Re 


Sn 


DMAA 


20 


20 


Pb 


La 


THF 


20 


10 


Pb 


W 


DMAA 


40 


20 


Cu 
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100 1256 
200 1056 
100 1069 
400 1110 
400 1082 

100 1068 

400 1851 

100 1078 

50 1118 

50 1252 

200 1040 

50 1088 

50 1086 

100 1073 

400 1926 

100 1085 

m X XT C I ii^ A 40 40 4oo iio6 

50 1201 

200 1460 
50 711 

lil Ni Sn THF 10 40 200 1272 

|'j S 1 ? P iG1 * 10 40 50 1059 

50 1718 

100 1060 

50 1292 

lil Ke Ce DMAA 5 40 50 1047 

Id £ c DMAA 20 10 200 792 

50 1057 

400 1045 

50 1074 

H " l r v umaa 20 5 400 1070 

100 758 

200 1047 

400 2009 

100 1695 

400 2255 

200 1701 

200 1366 

50 1073 

200 1090 

100 1677 

50 1421 

50 945 

400 453 

400 1303 

400 1102 
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[0034] 



The results in TABLE 4 are then subjected to an Analysis of Variance (ANOVA) 
analysis that includes the main effects and all the two-way interactions of the six 
factors (TM Type, TM Amount, LM type, LM amount, CS Amount, and CS Type). 
Results of the ANOVA are shown in TABLE 5. 
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The client and the investigator observe the rows of TABLE 5 that contain 
interactions. In the TABLE 5, only three of the interactions, marked **, show very 
strong evidence of statistical significance (P<0.001), and one, marked *, shows 
moderately strong evidence (P<0.02). Two show weak evidence (P-0.05). The rest 
show no evidence of interaction. The client and the investigator then adjust the 
weighted probabilities in the computer matrix according to the observed 
statistically significant results. The probabilities are increased for all the strong 
interactions and decreased for weak interactions. The following algorithm is used 
as illustrated in TABLE 6: (1) Very strong interaction: increase the matrix amount by 
half a distance to 1.0. (2) Moderately strong interaction: increase by .25 the 
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distance to 1.0. (3) Weak evidence: no change. (4) No evidence: decrease by half 
the distance to zero. 

TABLE 6 





I'M type 


TM 

Amount 


LMtype 


LM 

Amount 


CSType CS 

Amount 


TM type 


1 


0.8+.1 


0.3+.35 


0.3-.15 


0.3 


0.3-.15 


TM Amount 


0.8+.1 


1 


0.3+.35 


0.1-.05 


0.3-.15 


0.1-.05 


LM type 


0.3+.35 


0.3+.35 


1 


0.8 


0.3-.15 


0.1-.05 


LM Amount 


0.3-.15 


0.1-.05 


0.8 


1 


0.1-.05 


0.1+.225 


CSType 


0.3 


0.3-.15 


0.3-.15 


0.1-.05 


1 


0.8-0.4 


CS Amount 


0.3-.15 


0.1-.05 


0.1-.05 


0.1+.225 


,0.8-0.4 


1 



sg [0036] The revisions shown to TABLE 6, result in TABLE 7. 

|;g TABLE? 

%j TM type TM LM type LM CS Type CS 

*f. Amount Amount Amount 



S z 

n : s 


TMtype 


1 


.9 


.65 


.15 


0.3 


,15 


5 

Ii3 


TM Amount 


.9 


1 


.65 


.05 


.15 


.05 


3 ™ 
«^ 


LMtype 


.65 


.65 


1 


.8 


.15 


.05 




LM Amount 


.15 


.05 


.8 


1 


.05 


.325 


n 


CSType 


0.3 


.16 


.15 


.05 


1 


.4 




CS Amount 


.15 


.05 


.05 


.325 


.4 


1 



[0037] 



A full factorial experiment is organized and run according to the strongest 
interactions on the TM Type/TM Amount/LM Type variables (5x4x5 = 100 runs, 
fully replicated to 200 runs). Results are shown in TABLE 8. 
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TABLE 8 



TMType LMType CSType 



Fe 


V 


DMAA 


Fe 


w 


DMAA 


Fe 


Ce 


DMAA 


Fe 


La 


DMAA 


Fe 


Sn 


DMAA 


Cu 


V 


DMAA 


Cu 


w 


DMAA 


Cu 


Ce 


DMAA 


Cu 


La 


DMAA 


Cu 


Sn 


DMAA 


Ni 


V 


DMAA 


Ni 


W 


DMAA 


Ni 


Ce 


DMAA 


Ni 


La. 


DMAA 


Ni 


Sn 


DMAA 


Pb 


V 


DMAA 


Pb 


W 


DMAA 


Pb 


Ce 


DMAA 


Pb 


La 


DMAA 


Pb 


Sn 


DMAA 


Re 


V 


DMAA 


Re 


w 


DMAA 


Re 


Ce 


DMAA 


Re 


La 


DMAA 


Re 


Sn 


DMAA 


Fe 


V 


DMAA 


Fe 


W 


DMAA 


Fe 


Ce 


DMAA 


Fe 


La 


DMAA 


Fe 


Sn 


DMAA 


Cu 


V 


DMAA 


Cu 


W 


DMAA 


Cu 


Ce 


DMAA 


Cu 


La 


DMAA 


Cu 


Sn 


DMAA 


Ni 


V 


DMAA 


Ni 


W 


DMAA 


Ni 


Ce 


DMAA 


Ni 


La 


DMAA 


Ni 


Sn 


DMAA 


Pb 


V 


DMAA 





LM 


CS Amount TON 




t 


Amount 






5 


10 


100 


1 ns 

I I JO 


5 


10 


100 


1 1T7 


5 


10 


100 


1157 


5 


10 


100 


l*rZrf 


5 


10 


100 




5 


10 


100 




5 


10 


100 


Iftdft 


5 


10 


100 


1 1 50 


5 


10 


100 


1 17A 
1 1 /O 


5 


10 


1 w 




5 


10 


100 


QQA 


5 


10 


100 


fiOA 


5 


10 


100 


y\)D 


5 


10 


100 


K4fi 


5 


10 


100 


079 


5 


10 


100 




5 


10 


100 


0£5 


5 


10 


100 


JOD 


5 


10 


100 


700 


5 


10 


100 


190 


5 


10 


100 


549 


5 


10 


100 


767 


5 


10 


100 


491 


5 


10 


100 


726 


5 


10 


100 


511 


10 


10 


100 


1002 


10 


10 


100 


1038 


10 


10 


100 


1124 


10 


10 


100 


1211 


10 


10 


100 


1388 


10 


10 


100 


1000 


10 


10 


100 


1069 


10 


10 


100 


1064 


10 


10 


100 


1278 


10 


10 


100 


1269 


10 


10 


100 


1061 


10 


10 


100 


1136 


10 


10 


100 


977 


10 


10 


100 


1001 


10 


10 


100 


1487 


10 


10 


100 


1048 



Pb 
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10 
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10 
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10 
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10 
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10 
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10 
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Cu 


Sri 


DMAA 


40 


10 


100 


Ni 


V 


DMAA 


40 


10 


100 


Ni 


W 


DMAA 


40 


10 


100 



1188 
1333 
907 
1155 
1028 
839 
834 
1308 
1203 
879 
877 
888 
983 
759 

.0 wu v DMAA 20 10 100 1000 

" " 1016 

1146 
1236 
1205 
1149 
1062 
1289 
1374 
1668 

n « v dmaa 20 10 100 1126 

]t Pb W DMAA 20 10 100 1449 

!i " DU ^ ™' AA " 1476 

1592 
1828 

□ Ke v DMAA 20 10 100 1 136 

1728 
1481 
2336 
1928 
765 
741 
715 
475 
590 
1000 
1061 
1085 
1181 
1153 
1198 
1367 
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965 


1 A 
1U 


1 ft 

10 


100 


1480 


1 A 
1U 


10 


100 


1038 


1 A 
10 


10 


100 


1182 


1 ft 

10 


10 


100 


1016 


1 ft 

10 


10 


100 


979 


1 ft 

10 


10 


100 


828 


t ft 

10 


10 


100 


1204 


1 ft 

10 


10 


100 


1313 


20 


10 


100 


874 


20 


10 


100 


923 


20 


10 


100 


840 


20 


10 


100 


1017 


20 


10 


100 


700 


20 


10 


100 


1000 


20 


10 


100 


1046 


20 


10 


100 


1097 


20 


10 


100 


1172 


20 


10 


100 


1226 


20 


10 


100 


1106 


20 


10 


100 


1249 


20 


10 


100 


1201 


20 


10 


100 


1331 


20 


10 


100 


1302 


20 


10 


100 


1362 


20 


10 


100 


1308 


20 


10 


100 


1665 


20 


10 


100 


1558 


20 


10 


100 


1942 


OA 


1 ft' 

10 


100 


1390 


20 


10 


100 


1629 


20 


10 


100 


1731 


20 


10 


100 


2401 


20 


10 


100 


2327 


40 


10 


100 


748 


40 


10 


100 


674 


40 


10 


100 


714 


40 


10 


100 


691 
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S Ce £ 10 100 10 2« 

<-u Ce DMAA 40 io inn inio 

Cu La DMAA 40 10 J gg 

N W K 2 10 100 1258 

S Sn DMAA 2 '° ^ £ 

Ph v 40 10 "» 1663 

Pb V DMAA 40 10 ioo 1437 

Pb S„ DMAA 40 J 22 

Re V DMAA 40 10 JS JJJJ 

W DMAA 40 10 100 TOT 

Re Ce DMAA 40 10 00 2329 

Re La DMAA 40 10 22 ?S? 

Re Sn DMAA 40 10 22 ££ 



IB 



[0038] 



An ANOVA analysis of variance of the TABLE 8 data is illustrated in TABLE 9. 



Source 
TMType 



1 T%nw«MT. 10*51748 3483916 242.42 0 

^^^I yPC 16 1330652 133 «*2 83166 5 79 0 

U .^v£™ m0Un ! 12 22425009 22425009 1868751 130.03 2 

.H ^ ype ! TMAmount 12 1489975 1489975 124186 8 64 n 

1 - TMType *LMType*TMAmount 48 3450829 3450829 71892 55 2 

^ Error inn i^mia <^t<^ „ . u 



Total 





TABLE 9 




DF 


Seq SS 


Ad) SS 


4 


4777245 


4777246 


4 


2432949 


2432949 


3 


10451748 


10451748 


16 


1330652 


1330642 


12 


22425009 


22425009 


12 


1489975 


1489975 


48 


3450829 


3450829 


100 


1437139 


1437139 


199 


47795548 





Adj MS F p 



IMtw ^ 4777246 1194311 83.10 0 

!^JI ype 4 2432949 2432949 608237 4232 0 

™Amount 3 trusnaa iami-mq ,,0^, ? 



14371 



[0039] The ANOVA analysis detects a statistically significant 3-way interaction, which 
is a lead to high value formulations with high levels (TMA= 40) of Re in the 
presence of La or Sn. 

[0040] While preferred embodiments of the invention have been described, the 

present invention is capable of variation and modification and therefore should not 
be limited to the precise details of the Examples. The invention includes changes 
and alterations that fall within the purview of the following claims. 
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